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A pseudorandom point in an ergodic dynamical system over a computable metric space is a point 
which is computable but its dynamics has the same statistical behavior as a typical point of the 
system. 

It was proved in (2) that in a system whose dynamics is computable the ergodic averages of 
computable observables converge effectively. We give an alternative, simpler proof of this result. 

This implies that if also the invariant measure is computable then the pseudorandom points are a 
set which is dense (hence nonempty) on the support of the invariant measure. 



1 Introduction 

We will consider abstract algorithmic questions concerning the evolution of a dynamical system. In 
particular, the algorithmic estimation of the speed of convergence of ergodic averages and the recursive 
construction of points whose dynamics is typical for the system. 

This latter problem is related to the possibility of computer simulations, as actual computers can only 
calculate the evolution of computable initial conditions. 

Let X be a separable metric space and let the iterations of a map T : X — > X define a dynamics. Let p 
be an invariant measure for the dynamics (p(A) = p(T~ l {A)) for each measurable set A). The famous 
Birkhoff ergodic theorem says that if p is ergodic (the only sets which are invariant for the dynamics 
have full or null measure) then 

lim - Y\f(T n (x)) = [ f dp , p — almost everywhere. 

Similar results can be obtained for the convergence in the L 2 norm, and others. The above result tells 
that, if the system is ergodic, there is a fulll measure set of points for which the averages of the values of 
the observable / along its trajectory (time averages) coincides with the spatial average of the observable 
/. Such a point could be called typical for / . Points which are typical for each / which is continuous 
with compact support are called typical for the measure p (see Definition \Y2\. 

Estimating the speed of convergence. Many, more refined results are linked to the speed of convergence 
of the above limit. We consider this problem from the point of view of Computable Analysis. In the paper 
[2] some absttact results imply that in a computable ergodic dynamical system, the speed of convergence 
of such averages can be algorithmically estimated. On the other hand it is also shown that there are non 
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ergodic systems where this kind of estimations are not possible. In this paper we present a simple, direct 
way to prove the following result: 

Theorem A. If (X,ii,T) is an ergodic dynamical system and T is a.e. computable, then for each com- 
putable L 1 observable f the ergodic average A n {x) = ^L"=o f(T'(x)) converge effectively a.e. to ffdp, 
(see Definition \1 1 \for the precise definition of this kind of convergence). 

This theorem has some interesting consequences, as we are going to illustrate. 

Computable points having typical statistical behavior. The set of computable points in X (see Defi- 
nition [3) being countable, is a very small (invariant) set, compared to the whole space. For this reason, 
a computable point could be rarely be expected to be typical for the dynamics, as defined before. More 
precisely, the Birkhoff ergodic theorem and other theorems which hold for a full measure set, cannot help 
to decide if there exist a computable point which is typical for the dynamics. Nevertheless computable 
points are the only points we can use when we perform a simulation or some explicit computation on a 
computer. 

A number of theoretical questions arise naturally from all these facts. Due to the importance of 
forecasting and simulation of a dynamical system's behavior, these questions also have some practical 
motivation. 

Problem 1 Since simulations can only start with computable initial conditions, given some typical statis- 
tical behavior of a dynamical system, is there some computable initial condition realizing this behavior? 
how to choose such points? 

Such points could be called pseudorandom points. Meaningful simulations, showing typical behav- 
iors of the dynamics can be performed if computable, pseudorandom initial conditions exist ( and can be 
computed from the description of the system). 

We remark that it is widely believed that computer simulations produce correct statistical behavior. 
The evidence is mostly heuristic. Most arguments are based on the various "shadowing" results (see e.g. 
[9] chapter 18). In this kind of approach (different from ours), it is possible to prove that in a suitable 
system every pseudo-trajectory, as the ones which are obtained in simulations with some computation 
error, is close to a real trajectory of the system. However, even if we know that what we see in a 
simulation is near to some real trajectory, we do not know if this real trajectory is typical in some sense. 

The main limit of this approach is that shadowing results hold only in particular systems, having 
some uniform hyperbolicity, while many physically interesting systems are not like this. 

In our approach we consider real trajectories instead of "pseudo" ones and we ask if there is some 
computable point which behaves as a typical point of the space. Thanks to a kind of effective Borel- 
Cantelli lemma, in [6 ] the above problem is solved affirmatively for a class of systems which satisfies cer- 
tain technical assumptions which includes systems whose decay of correlation is faster that Clog 1 (time). 
In this paper we prove the following more general result, as a consequence of the above Theorem A: 

Theorem B. If (X,n,T) is a computable dynamical system as above and ii is a computable invariant 
ergodic measure, there exist computable points xfor which it holds: 

Um /M+/(n*)) + ...+/(r- 1 M) = I 

n-f*> n J 

for any continuous function f : X — > R with compact support. 
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The above theorem hence states that in such systems there are pseudorandom points. 
Physical measures and computability. To apply the above corollary to concrete systems the main 
difficulty is to verify that the invariant measure is computable. In [6] and (H it is shown that this is 
verified for the physicaQ invariant measure (the natural invariant measure to be considered in this cases) 
in several classes of interesting systems as uniformly hyperbolic systems, piecewise expanding maps and 
interval maps with an indifferent fixed point. On the other hand there are computable systems having no 
computable invariant measure ( [8 ]) this shows some subtlety in this kind of questions. 

The way we handle computability on continuous spaces is largely inspired by representation theory 
(see Ifl6l ,ir3l0. However, the main goal of that theory is to study, in general topological spaces, the way 
computability notions depend on the chosen representation. Since we focus only on computable metric 
spaces, we do not use representation theory in its general setting but instead present computability notions 
in a self-contained way, and hopefully accessible to non-specialists. 

2 Computability 

The starting point of recursion theory was to give a mathematical definition making precise the intuitive 
notions of algorithmic or effective procedure on symbolic objects. Several different formalizations have 
been independently proposed (by Church, Kleene, Turing, Post, Markov...) in the 30's, and have proved 
to be equivalent: they compute the same functions from N to N. This class of functions is now called the 
class of recursive functions. As an algorithm is allowed to run forever on an input, these functions may 
be partial, i.e. not defined everywhere. The domain of a recursive function is the set of inputs on which 
the algorithm eventually halts. A recursive function whose domain is N is said to be total. 

We now recall an important concept from recursion theory. A set E C N is said to be recursively 
enumerable (r.e.) if there is a (partial or total) recursive function (p : N — > N enumerating E, that is 
E = {(p(n) : n 6 N}. If E ^ 0, <p can be effectively converted into a total recursive function y/ which 
enumerates the same set E. 

2.1 Algorithms and uniform algorithms 

Strictly speaking, recursive functions only work on natural numbers, but this can be extended to the 
objects (thought as "finite" objects) of any countable set, once a numbering of its elements has been 
chosen. We will use the word algorithm instead of recursive function when the inputs or outputs are 
interpreted as finite objects. The operative power of algorithms on the objects of such a numbered set 
obviously depends on what can be effectively recovered from their numbers. 

1 In general, given (X,T) there could be infinitely many invariant probability measures. Among this class of measures, 
some of them are particularly important because they are related to what can be seen in real experiments. Suppose that we 
observe the behavior of the system (X,T) through a class of continuous functions /; : X — s- R. We are interested in the statistical 
behavior of /,- along typical orbits of the system. Let us suppose that the time average along the orbit of x exists 

A x (f i )=lim^f i (T»(x)). 

This is a real number for each /,-. Moreover A x (fj) is linear and continuous with respect to small changes of /,■ in the sup norm. 
Then the orbit of x acts as a measure fl x and A x ■(/,) = Jf, dfl x (moreover this measure is also invariant for T). This measure is 
physically interesting if it is given by a "large" set of initial conditions. This set will be called the basin of the measure. If X is 
a manifold, it is said that an invariant measure is physical if its basin has positive Lebesgue measure (see 1 17] for a survey and 
more precise definitions). 
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More precisely, let X and Y be numbered sets such that the numbering of X is injective (it is then a 
bijection between N and X). Then any recursive function q> : N — > N induces an algorithm s£ :X — > Y. 
The particular case X = N will be much used. 

For instance, the set Q of rational numbers can be injectively numbered Q = {qo,qi,...} in an 
effective way: the number i of a rational a/b can be computed from a and b, and vice versa. We fix such 
a numbering: from now and beyond the rational number with number i will be denoted by qt. 

Now, let us consider computability notions on the set K of real numbers, introduced by Turing in 
El. 

Definition 1 Let xbe a real number. We say that: 

• x is lower semi-computable if the set {/ € N : qt < x} is r.e. 

• x is upper semi-computable if the set {i € N : qt > x} is r.e. 

• x is computable if it is lower and upper semi-computable. 

Equivalently, a real number is computable if and only if there exists an algorithmic enumeration of a 
sequence of rational numbers converging exponentially fast to x. That is: 

Proposition 1 A real number is computable if there is an algorithm .sif : N — > Q such that \stf {n) — x\ < 

1~ n for all n. 

Uniformity. Algorithms can be used to define computability notions on many classes of mathematical 
objects. The precise definitions will be particular to each class of objects, but they will always follow the 
following scheme: 

An object O is computable if there is an algorithm 

$4 :X ->Y 

which computes O in some way. 

Each computability notion comes with a uniform version. Let (f?,)i S N be a sequence of computable 
objects: 

Oi is computable uniformly in i if there is an algorithm 

j#:NxX^F 
such that for all i, £/[:=£/ (i, •) : X — > Y computes Oj. 

For instance, the elements of a sequence of real numbers (xi)i e ^ are uniformly computable if there is 
a algorithm j</:NxN->Q such that \s/(i,n) -x t \ < 2~ n for all i,n. 

In each particular case, the computability notion may take a particular name: computable, recursive, 
effective, r.e., etc. so the term "computable" used above shall be replaced. 
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2.2 Computable metric spaces 

A computable metric space is a metric space with an additional structure allowing to interpret input 
and output of algorithms as points of the metric space. This is done in the following way: there is a 
dense subset (called ideal points) such that each point of the set is identified with a natural number. The 
choice of this set is compatible with the metric, in the sense that the distance between two such points 
is computable up to any precision by an algorithm getting the names of the points as input. Using these 
simple assumptions many constructions on metric spaces can be implemented by algorithms. 

Definition 2 A computable metric space (CMS) is a triple X = (X,d,S), where 

(i) (X,d) is a separable metric space. 

(ii) S = {si}ien is a dense, numbered, subset ofX called the set of ideal points. 

(Hi) The distances between ideal points d(s{,Sj) are all computable, uniformly in i,j (there is an algo- 
rithm s/ : N 3 ->• Q such that \s/(i,j,n) - d(si,Sj)\ < 2~ n ). 

S is a numbered set, and the information that can be recovered from the numbers of ideal points is 
their mutual distances. Without loss of generality, we will suppose the numbering of S to be injective: it 
can always be made injective in an effective way. 

We say that in a metric space (X,d), a sequence of points (x n ) ne ?$ converges recursively to a point x 
if there is an algorithm D : Q —> N such that d(x n ,x) < e for all n > D(e). 

Definition 3 A point x G X is said to be computable if there is an algorithm s/ : N — > S such that 
(s/ («)) n eN converges recursively to x. 

We define the set of ideal balls to be SB := {B(s^qf) : 57 G S,0 < qj G Q} where B(x,r) = {y G X : 
d(x,y) < r} is an open ball. We fix a numbering SB = {Bq,B\,. ..} which makes the number of a ball 
effectively computable from its center and radius and vice versa. SB is a countable basis of the topology. 

Definition 4 (Effective open sets) We say that an open set U is effective there is an algorithm s/ : 
N -)• SB such that U = \J n s/(n). 

Observe that an algorithm which diverges on each input n enumerates the empty set, which is then an 
effective open set. Sequences of uniformly effective open sets are naturally defined. Moreover, if (E/;),- e N 
is a sequence of uniformly effective open sets, then U ; [/; is an effective open set. 

Definition 5 (Effective Gg-set) An effective Gg-set is an intersection of a sequence of uniformly effec- 
tive open sets. 

Obviously, an intersection of uniformly effective Gg-sets is also an effective Gg-set. 

Let (X,Sx = {jf ,sf > —}idx) and (Y,Sy = {s\, s\ , be computable metric spaces. Let also Bf 

and Bj be enumerations of the ideal balls in X and Y. A computable function X — >■ Y is a function whose 
behavior can be computed by an algorithm up to any precision. For this it is sufficient that the pre-image 
of each ideal ball can be effectively enumerated by an algorithm. 

Definition 6 (Computable Functions) A function T : X — >■ Y is computable ifT~ 1 (Bj) is an effective 
open set, uniformly in i. That is, there is an algorithm s/ :NxN^ S3 X such that T ~ 1 (Bj ) = \J n s/ (/, n) 
for all i. 

A function T : X Y is computable on D C X if there are uniformly effective open sets Ui such that 
T- l (B Y i )C\D = U i C\D. 
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2.3 Computable measures 

Let us consider the space PM{X) of Borel probability measures over X. Let Co(X) be the set of real- 
valued bounded continuous functions onX. We recall the notion of weak convergence of measures: 

Definition 7 p n is said to be weakly convergent to p if J f dp n — > J f dp for each f G Cq(X). 

Let us introduce the Wasserstein-Kantorovich distance between measures. Let p\ and p% be two 
probability measures on X and consider: 



where l-Lip(X) is the space of functions on X having Lipschitz constant less than one. 
Proposition 2 (see [lj Prop 7.1.5) 

1. W\ is a distance and ifX is bounded, separable and complete, then PM{X) with this distance is a 
separable and complete metric space. 

2. IfX is bounded, a sequence is convergent for the W\ metrics if and only if it is convergent for the 
weak topology. 

Item (1) has an effective version: PM(X) inherits the computable metric structure of X. Indeed, 
given the set Sx of ideal points of X we can naturally define a set of ideal points S PM ( X ) i n PM(X) 
by considering finite rational convex combinations of the Dirac measures 8 S supported on ideal points 
s G Sx- This is a dense subset of PM(X). The proof of the following proposition can be found in ( fiXlD . 

Proposition 3 IfX bounded then (PM (X) ,W\ , Spm(x)) is a computable metric space. 

A measure p is then computable if there is a sequence p n G S PM (x) converging exponentially fast to 
p in the W\ metric (and hence p n weakly converge to p). 

2.4 Computable probability spaces 

To obtain computability results on dynamical systems, it seems obvious that some computability condi- 
tions must be required on the system. The "good" conditions, if any, are not obvious to specify. 

A computable function defined on the whole space is necessarily continuous. But a transformation 
or an observable need not be continuous at every point, as many interesting examples prove (piecewise- 
defined transformations, characteristic functions of measurable sets,... ), so the requirement of being 
computable everywhere is too strong. In a measure-theoretical setting, a natural weaker condition is to 
require the function to be computable on a set of full measure. It can be proved that such a function can 
be extended to a function which is computable on a full-measure effective Gg-set (see iflTlfTOlO . 

Definition 8 A computable probability space is a pair (X,p) where X is a computable metric space 
and p a computable Borel probability measure on X. 

Let Y be a computable metric space. A function (X,p) — » F is almost everywhere computable (a.e. 
computable for short) if it is computable on an effective Gg-set of measure one, denoted by dom/ and 
called the domain of computability of /. 

A morphism of computable probability spaces f : (X,p) — > (Y, v) is a morphism of probability 
spaces which is a.e. computable. 




sup 

/ei-Lip(x) 
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Remark 1 A sequence of functions f n is uniformly a.e. computable if the functions are uniformly com- 
putable on their respective domains, which are uniformly effective G$-sets. Observe that in this case, 
intersecting all the domains provides an effective G$-set on which all f n are computable. In the following 
we will apply this principle to the iterates f„ = T n of an a.e. computable function T : X —tX, which are 
uniformly a.e. computable. 

Remark 2 The space L l (X,jx) (resp. I?(X,ji)) can be made a computable metric space, choosing 
some dense set of bounded computable functions as ideal elements. We say that an integrable function 
f : X — > R is L 1 (X , /l) -computable if its equivalence class is a computable element of the computable 
metric space l}(X,\l). Of course, if f = g \l-a.e., then f is L i (X , \l) -computable if and only if g is. Basic 
operations on L (X,fl), such as addition, multiplication by a scalar, min, max etc. are computable. 
Moreover, if T : X — >■ X preserves }X and T is a.e. computable, then f — > foT (from L 1 to l}) is 
computable (see / Ii2 1/ ). 

2.4.1 Application to convergence of random variables 

Here, (X,ju) is a computable probability space, where X is complete. 
Definition 9 A random variable on is a measurable function f :X — >■ R. 

Definition 10 Random variables f n effectively converge in probability to f if for each e > 0, pt{x : 
\f n (x) — f(x)\ < e} converges effectively to 1, uniformly in £. That is, there is a computable function 
n(e,8) such that for all n > n(e,8), jJ.{\f n — f \ > e} < 8. 

Definition 11 Random variables f n effectively converge almost surely to f if f' n = sup )t> , 1 \fk—f\ ef- 
fectively converge in probability to 0. 

The following result (O, Theorem 2) shows that if a sequence /„ converges effectively a.s. to / then 
there are computable points which for which f„(x) — > f(x). 

Theorem 1 Let X be a complete metric space. Let f n ,f be uniformly a.e. computable random variables. 
If f n effectively converges almost surely to f then the set {x : f n (x) — > f(x)} contains an effective Borel- 
Cantelli set (see the Appendix for the precise definition). 

In particular, it contains a sequence of uniformly computable points which is dense in Supp(/l). 

Remark 3 Moreover, the effective Borel Cantelli Set found above depends algorithmically on f n and on 
the function n(S,s) giving the rate of convergence (see the proof of Theorem 2 in Hence the result 
is uniform in f„ and n(S,s). 

2.5 Effective L 1 , L 2 convergence 

Let (X,ju, T) be a computable measure-preserving system and / a L 1 -computable function (in the sense 
that it is a computable point of the metric space L 1 ). It was proved in [2 J that the (L 1 ,L 2 and almost sure) 
convergence of the Birkhoff averages of / is effective as soon as the norm of the limit f* is computable. 
Here we give an alternative proof in the ergodic case which is simpler as it uses the classical convergence 
result instead of giving a "constructive" proof. 

Let us call (X,fi,T) a computable ergodic system if (X,ju) is a computable probability space where 
T is an endomorphism (i.e. an a.e. computable measure-preserving transformation) and (X,fi,T) is 
ergodic. Let | \f\ \ denote the L 1 norm or the L 2 norm. 



14 



Speed of convergence of averages and pseudorandom points 



Proposition 4 Let (X,ii,T) be a computable ergodic system. Let f be a computable element ofL} (X,ju) 
(resp. L 2 (X,ii)). 

The L l convergence (resp. L 2 convergence) of the Birkhoff averages off is effective. 

Proof. Replacing / with f — f fdii, we can assume that / fd]X = 0. Let A„ = (f + / o T + . . . + 
f oT n ~ l )/n. The sequence ||A„|| is computable ( see Remark [2]) and converges to by the ergodic 
theorems. 

Given p e N, we write m£Nasm = np + k with < k < p. Then 



f n-\ 



A np+k = — !— - I £M p ° ? pi + kA k a 
np + K \ . =0 



\\An P +k\\ < — 

< \\A P \\ + ^ 



< \\A P \\ + 

it 

Let e > 0. We can compute some p = p(e) such that \\A p \ \ < e/2. Then we can compute some 
n(e) > § 11/11- The function m(e) := n(e)p(e) is computable and for all m > m(e), \ \A m \ \ < e. □ 

2.6 Effective almost sure convergence 

Now we use the above result to find a computable estimation for the a.s. speed of convergence. 

Theorem 2 Let (X,fi,T) be a computable ergodic system. If f is L 1 (X , \x)-computable, then the a.s. 
convergence is effective. 

This will be proved by the following 

Proposition 5 If f is L l (X , ^-computable, and \\f\\oc is bounded, then the almost-sure convergence is 
effective ( uniformly in f and a bound on \ \f\ 

To prove this we will use the Maximal ergodic theorem: 

Lemma 1 (Maximal ergodic theorem) For f G L x (X,ii) and 8 > 0, 

M({ropWf|>«})<i||/||i. 

n O 

The idea is simple: compute some p such that | \A P \ \ \ is small, apply the maximal ergodic theorem to 
g := A P , and then there is no, that can be computed, such that a{ is close to Af, for n > no- 
Proof. Let £,5 > 0. Compute p such that \ \Ap\ \ < 8e/2. Applying the maximal ergodic theorem to 
g :=A f p gives: 

At({sup|A«|>5/2})<£. (2.1) 

n 

Now, A 8 n is not far from a{: expanding Af, one can check that 

r U O T H U 

A " =A " + np ' 
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whereu = (p-l)f + (p-2)foT + ...+foTP- 2 . ||u|| .<2^||/|| 00 soifn>iio>4(p-l)||/||../5, 
then ||Af — a£||oo < 8/2. As aresult, if |a{(jc)| > 8 for some n > n , then |Af (jc)| > 8/2. From (12.11 ). we 
then derive 

At({sup \A{,\>8})<s. 

n>no 

As hq can be computed from 8 and e, we get the result. □ 

Remark 4 This result applies uniformly to a uniform sequence of computable L°°(X,n) observables f n . 

We now extend this to L 1 (X, n) -computable functions, using the density of L°°(X, ju) in L 1 (X, jj.). 
Proof, (of Theorem 12) Let s, 8 > 0. For M G N, let us consider f' M € L°° (X,jj.) defined as 

f'M-! min (/' M ) iff{*)>° 
Jm[X) \ max(/,-M) iff(x)<0. 

Compute M such that ||/ — TmIIi < Se. Applying Proposition \5\ to f' M gives some such that 

ft rff 

jU({sup„>„ \A J n M \ > 8}) < £. Applying Lemma[l]to f^ = f-f' M gives M({sup„ \A J n M \ > 8}) < £. As a 
result, Ai({sup„>„ |a{| > 25}) < 2e. □ 

Remark 5 Also Theorem\2\applies uniformly on an uniform sequence of computable L}(X,}i) observ- 
ables f n . 

Remark 6 We remark that a bounded a.e. computable function, as defined in Definition \8\ is a com- 
putable element ofL}{X,\x) (see l\12\l ). Conversely, if f is a computable element ofL 1 (X,jJ.) then there 
is a sequence of uniformly computable functions f n that effectively converge ]X-a.e. to f. 



3 Pseudorandom points and dynamical systems 

As said before the famous Birkhoff ergodic theorem says that in an ergodic system, the time average 
computed along \i -almost every orbit coincides with space average with respect to \i. 

If a point x satisfies Equation 11.11 for a certain /, then we say that x is typical with respect to the 
observable /. 

Definition 12 A point x is /I -typical if x is typical w.r.t. every continuous function f : X — > R with 
compact support. 

We remark that from now on we will suppose that X is a complete metric space. We will see that 
such jU —typical points exist in computable ergodic systems. First we give a result for L 1 observables. 

Theorem 3 If (X,n,T) is a computable ergodic system, f is L (X,ju) and a.e. computable then there is 
a uniform sequence Xj of computable points which is dense on the support of H such that for each i 

lim-£/(r"(x,))= If An. 

Proof. Apply fheorem[T]to the sequence of uniformly a.e. computable functions /„ = A J which con- 
verge effectively almost-surely by fheorem[2] We obtain that the set of points for which -£/(7 vl (.Jt B )) — > 
f /d/x contains a sequence of computable points, as in the statement. □ 

Let g be the point- wise limit of a sequence of uniformly computable functions /„ (defined jU-a.e., as 
in Remark©: Theorem[3]can also be proved to hold for the observable g, i.e. in a computable ergodic 
system there exists computable points for which the Birkhoff averages of g converge to fgd/j.. 

Since it is possible to construct a r.e. set of computable functions which is dense in the space of 
compactly supported continuous functions we can also obtain the following 
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Theorem 4 If (X,il,T) is a computable ergodic system then there is a uniform sequence x n of com- 
putable points which in dense on the support of il such that for each n, x n is il— typical. 

Proof. Let us introduce (following Q) a certain fixed, enumerated sequence of Lipschitz functions. 
Let J^o be the set of functions of the form: 

g s , r>e = \l-\d(x,s)-r\ + /e\ + (3.1) 

where s G S, r, e € Q and \a\ + = max{a,0}. 

gs AE is a Lipschitz functions whose value is 1 inside the ideal ball B{s,r), outside B(s,r + e) and 
with intermediate values in between. It is easy to see that the real- valued functions g Si ,r ,e k '■ X — > R are 
computable, uniformly in i,j,k. 

Let & be the smallest set of functions containing J^o and the constant 1, and closed under max, min 
and rational linear combinations. Clearly, this is also a uniform family of computable functions. We fix 
some enumeration vjf of & and we write g n for Vjr(n) G & . We remark that this set is dense in the set 
of continuous functions with compact support. 

By Remark [6] moreover these are computable elements of L l {X,p,), hence Theorem [2] applies uni- 
formly to these observables. This means that we can apply Theorem Q] uniformly on this sequence. By 
intersecting all effective Borel-Cantelli sets given by Theorem[[] since the intersection of a uniform fam- 
ily of effective BC sets contains an effective BC set (see Remark [3] and Proposition [6]) and such a set 
contains a sequence of computable points which are dense in the support of jJ, (see Theorem [5), which 
are typical for all observables in & (in the same way as in Theorem [3]). Since each continuous function 
/ with compact support can be approximated in the L°° norm by a function in & (and in particular the 
approximating function in & have values near the values of / at almost each point) the statement is 
proved. □ 

3.1 Conclusion and some open question 

We have seen that in computable ergodic systems, the speed of a.e. convergence of ergodic averages is 
computable, and if moreover also the invariant measure is computable then there are computable points 
which are typical for the statistical behavior of the system: the pseudorandom points. 

The assumption about the computability of the measure is not redundant with the assumption about 
the computability of the system because, as said before, there are computable systems having not com- 
putable invariant measures. It is also interesting to remark that there are systems where the map T is 
computable, for which we can consider a non computable ergodic measure having no pseudorandom 
points (an example can be constructed considering a Bernoulli shift where the symbols have non com- 
putable probability). 

Our results about the existence of these points, however, tell not much about the computational 
complexity which is necessary to find them. It would be of practical importance to have fast algoritms 
for this computation. 

All these questions are related to another general (and vague) question we like to cite, which is of 
great practical importance: why, many "naive" simulations of dynamical systems give reasonable results? 
or more precisely: under which assumptions the simulations of a dynamical system by a computer (which 
is a kind of discrete model for a continuous phenomena) give correct results? 

4 Appendix: effective BC sets 

We recall some results from [6] which are used in the proofs of the present paper. 
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Given a measurable space X endowed with a probability measure p., the well known Borel-Cantelli 
lemma states that if a sequence of sets A# is such that £jU(A*) < °° then the set of points which belong 
to finitely many A&'s has full measure. It holds that if the Ak are given in some "effective" way (and p is 
computable) then this full measure set contains computable points, which can be effectively constructed. 

Definition 13 A sequence of positive numbers a,- is effectively summable if the sequence of partial sums 
converges effectively: there is an algorithm stf : Q — > N such that if A{e) = n then Y*i>n a i — £ - 

For the sake of simplicity, we will focus on the complements U n of the A n . 

Definition 14 An effective Borel-Cantelli sequence is a sequence (U n ) n€ f$ of uniformly effective open 
sets such that the sequence p(X\U n ) is effectively summable. 
The corresponding effective Borel-Cantelli set is \J k f] n >k Un- 

Proposition 6 The intersection of any uniform family of effective Borel-Cantelli sets contains an effec- 
tive Borel-Cantelli set. 

Theorem 5 Let X be a complete Computable Metric Space and p a computable Borel probability mea- 
sure on X. 

Every effective Borel-Cantelli set R, contains a sequence of uniformly computable points which is 
dense in the support of p. 
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